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1 Introduction 



The origin of the baryon asymmetry of the universe (BAU) is still an open question 
in cosmology and particle physics. New measurements of the cosmic microwave 
background, combined with large-scale structure data, yield a baryon to entropy 
ratio of PP 

r^B = ^ = (8.7±0.3) X 10-1^ (1) 

Three necessary conditions, stated by Sakharov j2], have to be fulfilled for a dy- 
namical generation of the baryon asymmetry: baryon number violation, C and CP 
violation, and departure from thermal equilibrium. In principle the standard model 
(SM) contains all these requirements, and the electroweak phase transition (EWPT) 
provides a natural mechanism for baryogenesis The baryon asymmetry is gen- 
erated during the phase transition by electroweak sphaleron processes. To avoid 
subsequent baryon number washout, the sphaleron rate has to be suppressed after 
the phase transition. Hence the transition must be strongly first order, i.e. the ex- 
pectation value of the Higgs field must be larger than about the critical temperature. 
In the SM there is no first-order phase transition for Higgs masses larger than about 
80 GeV 13], far below the current experimental bound of ttih > 114 GeV The 
SM therefore fails to explain the baryon asymmetry. Moreover, the CP violation in 
the CKM matrix is too small to produce a sufficiently large baryon number [U]. 

Over the years there have been many proposals to realize electroweak baryogen- 
esis in extended models (see, for instance, ref. ^ for a review). In supersymmetric 
theories, for example, a strong first-order PT can occur if the superpartner of the 
top quark is lighter than about 175 GeV [8J, and the baryon asymmetry can be 
generated by chargino transport jHlEj- Alternatively, the phase transition can be 
strengthened by the presence of SM singlets in the Higgs sector ^T]. A more gen- 
eral effective field theory approach can also be followed; there the Higgs sector is 
augmented by dimension-six operators to induce a first order phase transition and 
to provide additional CP violation [T21 CHI HI] • 

In this paper we revisit electroweak baryogenesis in the two-Higgs doublet model 
(2HDM) [lailiinilliliniliniin]- in addition to the SM Higgs, the 2HDM 
contains two extra neutral and charged Higgs particles. If these extra states couple 
sufficiently strongly, their thermal loop corrections can induce a strong first-order 
phase transition I17| IT^ I21j . In addition, a complex mass term, mixing the two 
Higgs doublets, provides a new source of CP violation, which can fuel baryogenesis. 

We examine the EWPT in the 2HDM, including explicit CP violation, using the 
finite temperature effective potential at one-loop order. In agreement with ref. [23], 
we find a strong phase transition for light Higgs masses of up to at least 200 GeV. 
The extra Higgses have to be heavier than about 300 GeV, depending somewhat on 
the model parameters. Turning on the CP-violating phase makes the phase transi- 
tion slightly weaker. We determine the profile of the bubble wall, which separates 
the broken and symmetric phase. Except for the case of very strong phase transi- 
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tions, we typically find thick bubble walls. The bubble wall is characterized by a 
varying complex phase between the two Higgs vevs. CP-violating interactions of the 
particles in the hot plasma, in particular the top quark, with the phase boundary 
then lead to different semiclassical forces acting on particles and antiparticles. Since 
we are dealing with thick bubble walls, we can apply the standard WKB formal- 
ism to compute the CP-violating source terms that enter the transport equations 
of electroweak baryogenesis Here we use the formalism recently laid out in 

ref. ^3], which makes sure that the correct dispersion relations of the Schwinger- 
Keldysh formalism ^S] are reproduced. Also a finite VT-scattering rate is included 
in the transport equations, which previously was set to equilibrium. 

We find that a wide parameter range allows for the generation the observed 
baryon asymmetry. Since the model contains only a single CP-violating phase, 
we can predict the electric dipole moments of the neutron and electron. They 
are typically found to be below the current experimental bounds, but should be 
detectable in next-generation experiments. 

2 The effective Higgs potential 

In its most general form, the 2HDM suffers from flavor changing neutral currents 
at tree- level. To avoid this, a discrete symmetry. Hi — > —Hi, =F(if (the other 

fields do not transform) , is usually invoked [24J , making sure that at most one Higgs 
doublet couples to the up- and down-type quarks, respectively. In the "— " case 
("type H"), the down-type quarks couple only to Hi, while the up-type ones couple 
to H2. In the other case ("type I"), Hi does not couple to the fermions at all. In the 
following only the coupling of the top quark will be relevant, so that we do not need 
to actually distinguish between types I and II. In section 4, where we will discuss 
electric dipole moments, we will focus on the type II case. 
The most general potential is [221 

Vo{Hi,H2) = ~filHlHi-i2lHlH2-fil{e"^HlH2 + h.c.) 

+^{HlHif + ^{HlH.f + \,{HlHi){HlH,) 

+X,\HlH,\' + ^ [{HlH,)' + h.c.) (2) 

and the Yukawa interactions read 

Cy = yH^Qzt' + h.c. + . . . (3) 

Without loss of generality the couplings Aj and the mass parameters /Xj can taken to 
be real. The mass term ^\e^^ breaks the aforementioned Z2 symmetry softly, without 
reintroducing tree-level flavor violation j^H]- It can be complex, in which case the 
Higgs potential breaks CP. In total, the Higgs potential contains 9 parameters, which 
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are 3 squared masses, 5 couplings, and 1 phase. One parameter is fixed by the Z- 
boson mass, leaving an 8-dimensional parameter space. The potential has to be 
bounded from below, which at tree-level translates into the constraints ^HI 

Ai > 0, Aa > 0, VA1A2 + A3 > 0, VA1A2 + A3 + A4 ± As > 0. (4) 

Let us first consider the case </> = 0, i.e. the Higgs potential without CP violation. 
It has been shown that in this case no charge breaking minima exist, provided 
the charged Higgs mass squared is positive Later on we will assume that 

this result generalizes to the one-loop level, including a small CP-violating phase. 
We can therefore restrict ourselves to the neutral fields, which we parameterize as 
Reif° = h\ and R.eH2 = h2- The potential then reads 

Voih, h^) = -fijhl - fijhl - 2filhh2 + + ^hl + (A3 + A4 + h)hlhl. (5) 
In the following we focus on the somewhat simpler case 

/ii = /i2, Ai = A2. (6) 

Moreover, this choice is favorable to generate large Higgs expectation values in the 
broken phase ^Hj. The Yukawa interaction (jH)) does not preserve these relations at 
the loop-level. At tree-level, eq. © implies the symmetry 

Hi ^ Hi (7) 

so that the minimum is at tan (3 = (/la) / (hi) = 1. With (hi) = (/12) = h = 123 GeV 
the extremal condition is then given by 

-fij + (Ai + A3 + A4 + As)/^^ = 0. (8) 

The mass matrix is block-diagonal and we obtain, besides 3 massless Goldstone 
bosons, 5 physical Higgs bosons. They consist of a pair of charged Higgses 
neutral scalars and H^, and a pseudoscalar A^, with the corresponding squared 
masses as follows: 



m^± = 2/^2 - 2(A4 + A5)/l^ (9) 

mlo = 2i4-A\5h'^, (10) 

m^o = 2^2 _ 2(_A^ + A3 + A4 + A5)/^^ (11) 

mlo = 2(Ai + A3 + A4 + As)/^'. (12) 



These relations can be used to define the model in terms of /x| and the 4 Higgs 
masses. 

study of symmetry breaking in a general 2HDM has recently been presented in ref. j28| . 
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In the case of non-vanishing 0, CP is broken. We now parametrize the neutral 
Higgs fields as 

Hi = hre-'^\ Hi = h2e'^\ (13) 

Note that the potential only depends on the combination 9 = 9i + 92- In the 
minimum we can always choose the gauge such that 9\ = 92 = 9 j2. Still assuming 
the relations the potential of the neutral fields reads 

^0(/il,/i2,^) = -/i?(/l? + /l2) -2/i^/ll/l2COs(^ + 0) + y(/lt + /l^) 

+(A3 + A4 + A5Cos(2^))/i?/i^. (14) 

Using the notation (0) = ^9, we obtain two extremal conditions, 

- /ii - /i3Cos(^9 + 0) + (Ai + A3 + A4 + A5Cos(2'i9))/i^ = 

/i^sin(^9 + 0) - A5sin(2?9)/i2 = 0. (15) 

The squared Higgs boson masses take the form 

m\± = -2/i? + 2(Ai + A3)/^^ 

= -/i? + 2(Ai + A3 + X4)h^ + \/f^f + 4A5 cos(2^?)/if/i2 + 4Xlh\ 
m]j^ = -2/i2 + 4Al/l^ 

= -1^1 + 2(^1 + ^3 + A4)/i' + 4A5 cos{2d)ii\h^ + A\lh\ (16) 

Note that the neutral Higgs states are now mixtures, with scalar and pseudoscalar 
contents. Again, these relations can be inverted to parameterize the model in terms 
of the Higgs masses, ^\ and 0. 

At zero temperature the one-loop contribution to the effective potential is given 

by 

V.{H,^H2) = Y.±^,mtlnf,. (17) 

i 

where m? = mf{Hi, H2) are field dependent mass eigenvalues, is the correspond- 
ing number of degrees of freedom, and " + (—)" applies to bosonic (fermionic) contri- 
butions, respectively. We choose Q = 246/ \/2 GeV for the renormalization scale. We 
take only the heaviest bosons, i.e. rrii = mH±,'mH2^^Hi {P'H± = = = 1), 
and the fermion with the largest Yukawa coupling, i.e. rrii = rrit [rit = 12), into 
account. For the top quark mass we have mf = y^^i^l-^s- All other particles can be 
safely neglected, owing to their small contributions to the effective one-loop poten- 
tial. 

We add counter-terms to the potential, such that the tree-level minimum and 



4 



Higgs masses are preserved at the one-loop level. This can be achieved by 



Vct{HuH2) = -6fil{HlH, + HlH2)-6fil{e''f'HlH2 + h.c.) 

+5\,\H\H^\' + ^ {{HlH.f + h.c.) . (1^ 



As already mentioned, the symmetry (|7j) no longer holds at the one-loop level, which 
we take care of by using 5\i 7^ 8X2- Three renormalization conditions are evidently 
given by 



diV^ + V^ct) 



d{Vi + Vct) 



dho 



d{Vi + Vct) 



89 



h-^=h2—h 



0, (19) 



meaning that the minimum of the potential V = Vq -\-Vi -\- Vct does not change 
with respect to the tree-level case. Preserving the values of the Higgs masses, which 
we compute from the second derivatives of V , provides another four conditions. 
So the coefficients of Vct are fixed. Since the conditions related to Higgs masses 
include non-linearities, the resulting equations for the counter-terms have to be 
solved numerically. 

The 2HDM is subject to a number of experimental constraints. In the consid- 
ered parameter range, the lightest Higgs boson is SM-like, and therefore has to obey 
the LEP bound of niH > 114 GeV 5J. The 2HDM does not respect the custodial 
symmetry of the SM. So there is the danger of large corrections to the electroweak 
precision observables. These corrections can be approximately described in terms of 
contributions to the self-energies, the so called "oblique" corrections. The relevant 
expressions for the 2HDM with CP- violation can be found in ref. To be consis- 
tent with observations, the mass splittings between the extra Higgs states should not 
be much larger than the W-mass. Later on we will set these masses equal to reduce 
the number of parameters. Oblique corrections then are automatically small. An- 
other important constraint comes from b — >■ S7, which in the type II model requires 
mH± ^ 200 GeV jHOl^- Constraints from the muon anomalous magnetic moment 
and from tau decays [32] are not relevant for the low values of tan (3, which we 
consider. 

At finite temperature the one-loop contribution to the effective potential is given 

by 

Vl = T'Y, us fs (^) + T' J2 u, f, (!^) , (20) 

B F 

where nB{F) counts the positive degrees of freedom for bosons (fermions). In the 
^One can also consider constraints from B^^ — mixing 
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high temperature hmit, m/T <^ 1, one obtains 



m\ —TT^ im? Trfi , / iv? \ 



\t) 90 ^ 24r2 127rT3 QAn^T^^"^ \cbT^ J ^^^^ 

'vry 720 +48T2 + 647r2T4'''l^Ci.r2y'' ^^^^ 

with Ci? = 7r^exp(3/2 — 27e) ~ 13.94 and cb = IQcp, and in the low temperature 
hmit, when m/T is large, 

In this low temperature limit the contributions from bosons and fermions have the 
same asymptotic behavior. 

We use these approximations because they are much more convenient to handle 
than the full integral expressions of ref. It turns out, however, that these 

limiting cases are not sufficient since some states cross from the high temperature to 
the low temperature regime. For an expression to be valid in the whole temperature 
range, we therefore use a smooth interpolation between the low- and high-T limits. 
For bosons we use eq. ()2H1 for m/T < 1.8 and eq. ()23|1 for m/T > 4.5, and for 
fermions eq. (j^^ for m/T < 1.1 and eq. (j^^ for m/T > 3.4. The interpolations 
are made in such a way that the functions as well as their derivatives match at the 
connecting points. The deviation between our approximation and the exact solution 
is less than 4%. Finally, the effective potential is given by 

Kff = + 1^1 + l^cT + V^. (24) 

In Vy we also take into account the contributions of l^-bosons [nw = 6) and Z- 
bosons {nz = 3). In perturbation theory the strength of a strong phase transition 
would be underestimated by resummation of the gauge boson contributions [33]. 
Therefore we do not resum these corrections to compensate for non-perturbative 
effects. In ref. j2l] the thermal contributions of the heavy Higgs bosons have been 
resummed, using the high temperature approximation for the thermal Higgs masses. 
We find that in the broken minimum the high temperature approximation is often 
not justified for the heavy Higgs states. Ignoring this fact, and nevertheless resum- 
ming the Higgs contributions by using thermal masses instead of the bare masses in 
eq. ()20p . we find a phase transition less than 15% weaker than in unresummed case. 
The effect is marginal for large heavy Higgs masses and becomes stronger for smaller 
ones. Including the light Higgs and the Goldstone bosons has an even smaller effect. 
In the results we present below, the Higgs contributions are not resummed. 



3 The phase transition 

The dynamics of the EWPT is governed by Ves{hi, h2,6,T). The critical tempera- 
ture, Tc, of a first-order phase transition is defined by the condition that the effective 
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potential has two degenerate minima, the symmetric minimum at (/ii)t = (^2)t = 
and the broken minimum at {hi)^ = fi > and (/i2)t = f2 > 0, which are 
separated by an energy barrier. The total Higgs expectation value we define by 
Vc = \f^\/ + where the factor \f2 is due to our normalization of the Higgs 
fields. Somewhat below Tc, at the nucleation temperature T„, bubbles of the broken 
phase start to nucleate and expand. Baryogenesis takes place outside the bubbles 
in the symmetric phase while, inside the bubbles, the sphaleron rate that provides 
(B + L)-violating processes has to be switched off. Otherwise the baryon asymme- 
try will be washed out after the phase transition. In order to preserve the created 
baryon asymmetry, the washout criterion 

e=|>l (25) 

must hold, i.e. the phase transition has to be sufficiently strong. 

In the following we analyze the parameter space with respect to the strength ^ 
of the phase transition. We focus on the case of degenerate heavy Higgs masses, 
which reduces the dimension of the parameter space. As input parameters we take 
/X3, 0, m/j = m/^j, and mn = niu^ = itlh^ = mH±. One finds that for larger values of 
(0 = 0.4 for example) the first-order phase transition can change into a two-stage 
transition if the heavy Higgs mass is sufficiently small. The potential then shows an 
additional local minimum. The phase transition proceeds by a second-order phase 
transition from the symmetric phase to this extra minimum, followed by a first- 
order phase transition to the low temperature broken phase®. We exclude these 
values from the parameter space and only define the strength ^ in the case of a pure 
first-order PT. 

Another important property that enters the transport equations discussed in 
chapter 5 is the wall profile of the expanding bubbles. If the nucleating bubbles 
have reached a sizable extent and expand with constant velocity, we can boost into 
the rest frame of the bubble wall and assume a planar wall. In principle one has to 
numerically solve the field equations of the Higgs fields, using an algorithm such as 
the one recently proposed in ref . ^7 • To achieve a sufficiently strong phase transition 
we are led to rn^j ^ m\. The effective potential is then characterized by a valley, 
corresponding to the single light Higgs field. During the phase transition the fields 
will follow this valley very closely, in order not to feel the heavy Higgs masses. So 
we can approximate the phase transition by single field dynamics. Numerically we 
determine the valley by minimizing the thermal potential at T^. with respect to /t2 
and 9 at fixed values of hi between the symmetric and broken phase. For a simple 
<y9^ model, with one real scalar field (p and a broken minimum at fc, the wall profile 

^The second-transition is in general too weak to avoid baryon number washout. Also, the non- 
zero Higgs vev outside the bubbles would lead to a suppression of the sphaleron rate and therefore 
also of the baryon asymmetry. So the two-stage transition does not allow us to generate the baryon 
asymmetry. 
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is exactly described by a kink solution, 

cpiz) = I (l - tanh-^) , (26) 

with a wall thickness L„ = ^^/vfJ{8Vb), where H is the height of the potential 
barrier and z is the coordinate orthogonal to the wall. We use this approximation 
for Lw, determining Vb as the maximal height of the potential barrier along the 
valley connecting the two minima. 

Let us now briefly discuss the behavior of and ^ with the input parameters. 
We require > 115 GeV to be consistent with the LEP bound on the Higgs 
mass. For increasing rriH and keeping the other parameters fixed, the wall thickness 
decreases, while the PT becomes stronger. This somewhat counter-intuitive result is 
due to the fact that the larger Higgs masses (fTBj) come from larger quartic couplings. 
So this limit actually does not lead to the decouphng of the heavy states. At some 
point perturbation theory will finally break down. Later on, when computing the 
baryon asymmetry, we will face another constraint. The gradient expansions is 
justified only for thick bubble walls, so we will require L^Tc > 2. In practice this 
leads to an upper bound on mn similar to the perturbativity constraint. 

In figures HHSl constant lines of ^ and are shown in the dependence of rrih and 
mn for different values of /i| and 0. The influence of the CP-violating phase on ^ 
and L„ is rather small, as can be seen from figs. ^ and El For small values of ^, or 
large values of L^, the lines are marginally shifted upwards. This behavior continues 
for increasing 0. If we choose = 0.4 the effect is negligible above ^ ~ 1.5, but 
below C, ~ 1.3 the PT changes into a two-stage one. The effect of increasing /i| is a 
shift to higher values for ttih. The comparison of /x| = 10000 GeV^ and 20000 GeV^, 
e.g. figs. El and El shows that the range of niH is moved to higher values, while the 
extent shrinks by around 20 GeV. Using a larger value of fi^ means that the same 
quartic couplings lead to heavier Higgses. The strength of the phase transition is 
more governed by the size of the quartic couplings than by the actual value of mn- 
Notice that here the bound on the charged Higgs mass from b — >■ 57, which we 
discussed above, is automatically satisfied in the case of a strong phase transition. 

In figs, m and El the one-loop corrections A = max|(5Aj/Aj| to the quartic cou- 
plings, i.e. the size of the counter-terms with respect to the tree-level terms, range 
from 15% for mn = 300 GeV to 50% for mn = 440 GeV. This means that in the 
case of a very strong phase transition, perturbation theory starts to break down 
and sizable corrections to our results have to be expected. For ^ ~ 1 higher-order 
corrections are well under control. These results agree with the findings of ref . [21] . 
In conclusion, we find that a wide range of parameters fulfills the requirements of 
electroweak baryogenesis. 

Up to now we have only discussed observables involving the fields hi and /i2- 
However, also the CP-violating phase 6, which varies along the bubble wall from 
6'sym to ^brk, is essential for baryogenesis. According to the above discussion, we 
compute the ^-profile approximately by minimizing the thermal potential at with 
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Figure 1: Lines of constant ^ and L„ in the m/j-m/^-plane for /i| = 10000 GeV^ 
and = 0. In addition, the hne of the relative size of the one- loop corrections 
A = max |5Aj/A,;| = 0.5 is shown. The Higgs masses are given in units of GeV. 
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Figure 2: The same plot as in fig. ^ but for the set /ig = 10000 GeV^ and = 0.2. 



9 



360 



120 



130 



140 



150 



160 



170 



180 



190 



Figure 3: The same plot as in figs. Hand but for fxj = 20000 GeV^ and = 0.2 

respect to /i2 and 6 at fixed values of hi between the symmetric and broken phase. 
As a representative example we show in fig. the ^-profile parametrized by hi for 
the set fil = 10000 GeV^ and = 0.2. There, 6 changes from ^sym = -0.29 to 
^brk = —0.06, which is indicated by the dotted lines. In a simplified manner, we 
describe the ^-profile by a kink ansatz, i.e. 



using the derived wall thickness and A6 = ^^brk — ^sym- The CP violation in the 
Higgs sector gives rise to complex fermion masses, which change while the particles 
pass through the bubble wall. We only take into account the top quark as the 
heaviest fermion. With our parametrization ()13|1 of the neutral Higgs components, 
one finds for the complex top mass 



where Pt is the angle between hi and h2 at T^, i.e. tan/^T = V2/V1, which is less, 
but rather close to 1. The top Yukawa coupling yt is chosen such that the top 
mass at zero temperature is 173 GeV. In general, the change in ^2(= dt) along the 
bubble wall is given by A92 = A6'/(l + tan^/^r), assuming that tan/?T is constant 
along the wall [38j. So there is an additional suppression for A62 for large tan/3. 
A4t enters the computation of the baryon asymmetry; in particular, the derivative 
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Figure 4: The phase 6 and the difference for the set = 10000 GeV^ = 0.2. 

(a) The change of 6 during the PT, as a function of hi (given in GeV), at fixed 
TUh = 150 GeV, tuh = 350 GeV. 

(b) A9 versus rrih (given in GeV) for mn = 330 GeV (solid) and 400 GeV (dashed). 



of Ot{z) induces the CP-violating the source term. Therefore a large value of 
enhances the baryon asymmetry. As shown in fig. EJd, A^ strongly depends on m^. 
Raising ruh or mn does increase the change in 6. We also find that A^ depends 
almost linearly on the coupling 0, whereas the influence of fi^ is small. Notice that 
A^ can be larger than the input phase 0. 



4 Electric dipole moments 

CP violation induces electric dipole moments (EDMs). The latest experimental 
limits for the neutron jSHl and electron [101 EDMs at 90% confidence level are 

\dn\ < 3.0 X 10"^^ e cm, (29) 
141 < 1.6 X 10"^^ e cm. (30) 

In the standard model the only source of CP violation originates from the Kobayashi- 
Maskawa matrix in the quark sector. Contributions to the EDMs arise first at the 
three-loop level, which results in a natural suppression, several orders of magnitude 
below current bounds. EDMs are therefore an ideal probe of new physics. 

In the 2HDM, EDMs are induced by scalar-pseudoscalar mixing in the neutral 
Higgs sector. The contributions to the EDMs can be computed in terms of param- 
eters Im(Z), which measure the degree of CP non-conservation and which are the 
imaginary parts of Higgs fields normalization constants jHj. The four CP- violating 
parameters 

Im(Zoi), Im(Zoi), lm{Zu), lm{Z2i), (31) 
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where i indicates each of the four neutral Higgs bosons, enter the calculation of the 
EDMs. They can be expressed in terms of components of the neutral Higgs mass 
matrix eigenvectors. The Goldstone boson does not contribute to these factors, since 
the corresponding Z's are real. Thus, the sum can be restricted to the three massive 
neutral bosons. Note that the parameters respect in addition the sum rules jH] 

J2 MZo^) = Yl I^(^o*) = MZu) = J2 = 0, (32) 

i i i i 

which means that CP violation vanishes if the masses of the neutral Higgs bosons 
are degenerate. 

In the 2HDM the dominant contributions to the electron EDM are two-loop 
amplitudes, which were first computed by Barr and Zee |42j. They demonstrated 
that the effect is enhanced with respect to the standard one-loop contributions 
Further two- loop diagrams, including the W-boson, were taken into account 
in investigations by Gunion and Vega ^Ij, Chang et al. [IS], as well as Leigh et 
al. jinj. In this work we use the results of Chang et al., ignoring some minor 
corrections discussed in ref. llHj- We end up with the following contributions 



dje = (4/e)fTop + (4/e)Sp + (4/e)S-lVp 

+ (4/ e)w-]oop + {del e)S2Lop + (4/ e)Sf loop- (33) 

When computing the EDM of the neutron one has to deal with hadronic effects, 
which make its relation to the partonic EDMs difficult. Various proposals have been 
made in the literature how to perform this calculation (see ref. jU] for a recent 
review). The dominant contributions to the neutron EDM come from the color 
EDMs (CEDMs) of the constituent quarks 4, k = u,d 



i ~ - 1 - _ ^ 



2 " " ^ 2 
and from Weinberg's three-gluon operator |3T] 



£ D iw/'^''"G;^^G"^^'^G'^'". (35) 
3 



The QCD-corrected coefficients du, 4 and w are given by 2- loop calculations 
l^m 1^ . Using the results of |17], based on QCD sum rule techniques, the neutron 
EDM reads 

(4/e)(4, 4) = (1 ± 0.5) (0.554 + 1.14) (36) 

and 

|(4/e)(w)| = 22 MeV (37) 

Thus there is an error of about 50% in (4/e)(4,4) and furthermore an error of 
about 100% in (4/e)(w). Moreover, it is not possible to determine the sign of 
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Figure 5: Lines of constant neutron (solid) and electron EDMs (dashed) for the 
set i^l = 10000 GeV^, (p = 0.2, dn is given in units of 10~^^ ecm, in units of 
10~^^ ecm, and Higgs masses in GeV. The lower dotted line indicates the bound 
^ = 1, the upper one — 2. 



{dn/e){w). Fortunately, this latter contribution turns out to be small (typically a 
1% correction). The same is true for the contributions of the quark EDMs du and 
dd (typically a 10% correction). 

Let us now discuss the relevance of the electron and neutron EDMs for the 
2HDM. One finds that in the analyzed parameter region the value of de is about 
five to thirty times smaller than the experimental limit of 1.6 x 10~^^ ecm. Thus, 
there emerges no additional constraint on the parameters. Let us focus on the 
importance of the different contributions to d^ and on the dependence of d^ on 
the input parameters. Since an EDM arises because of CP violation, we expect a 
larger value for an increasing CP phase 0. Indeed we find that de approximately 
doubles if we change from 0.2 to 0.4. Also raising yu| enhances de- Concerning the 
single contributions to de the largest ones originate from the top- and W^-loops, with 
(de) VK-ioop > whereas ((ie)t-ioop < 0- The absolute value of ((ie)t-ioop is somewhat 
smaller, but of similar magnitude as {de)w -loop- So the sum is a factor of about 
5-10 smaller than each individual contribution, and is then of the same order of 
magnitude as the Goldstone-loop contribution. Thus, all three parts are important 
for the electron EDM. We observe this behavior in the whole analyzed parameter 
region. We also investigate the dependence of de on the Higgs masses. For increasing 
both, rrih and uih, the value of de decreases. This tendency becomes apparent in 
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fig. Elwhere we compare lines of constant electron and neutron EDMs in the rrih-mH 
plane. We find that the larger the weaker is the dependence on the heavy Higgs 
mass. 

Similar to the electron EDM, the one of the neutron also lies below its experimen- 
tal bound of 3.0 x 10~^^ e cm in the analyzed parameter region. But in contrast to 
de it almost reaches this experimental limit. However, note that dn has quite a large 
error of about 50%, as pointed out. The limit of 3.0 x 10~^^ ecm is just saturated 
in the case of /Xg = 20000 GeV^ and (p = 0.4, for small Higgs masses. However, since 
the error band is large, there actually arises no constraint. For larger values of 
or 0, the neutron EDM of course increases and may exceed the measured bound in 
a wider mass range. The dependence of the neutron EDM on the input parameters 
is quite similar to that of the electron EDM. The lines of constant dn run approxi- 
mately parallel to those of de in the rrih-mH plane; the slope is just a little fiatter. 
We also find roughly a doubling for a change in from 0.2 to 0.4. The dominant 
contribution arises from the color EDM of the down-quark, which is about a factor 
3.5 larger than the one due to the up-quark CEDM. The part |(i„(w)| arising from 
the three-gluon operator is roughly an order 1% correction and can therefore be 
neglected. In summary, for the considered parameter ranges, both the electron and 
neutron electric dipole moments lie below the experimental limits. The value of de 
is about one order of magnitude below the observational bound, and because of the 
large error in the theoretical determination of dn, it also does not definitely exceed 
the bound set by experiments. 

5 Transport equations 

In this section we discuss the evolution of the particle distributions during the phase 
transition. The CP-violating interactions of particles in the plasma with the bubble 
wall create an excess of left-handed quarks over the corresponding antiquarks. This 
excess diffuses into the symmetric phase, where the left-handed quark density biases 
the sphaleron transitions to generate a net baryon asymmetry. 

Using the semiclassical WKB formalism El 122]! we obtain different dispersion 
relations for particles and antiparticles in the space-time dependent background of 
the Higgs expectation values. The dispersion relations then lead to force terms in the 
transport equations. The WKB method is justified when the de Broglie wavelength 
of the particles in the plasma is much shorter than the bubble wall thickness |22] ■ 
Hence the condition L^T ^ 1 has to be satisfied to legitimate an expansion in 
derivatives of the background Higgs fields. As demonstrated in section 01 we find 
that a large part of the parameter space does fulfill this condition. 

In the 2HDM, baryogenesis is driven by top transport. So we can focus the 
discussion on the case of a single Dirac fermion, with a space-time dependent mass 
ReA^(z) + i'j^lmAi{z), where Ai{z) = m{z)e^^^^\ The dispersion relation to first 
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order in gradients is given by [TU 



E = Eo±AE = EoT s^^^, (38) 



where E^ = a/p^ + and £"02 = ^/j^l-rn^ in terms of the kinetic momentum. 
The prime denotes the derivative with respect to z, and the upper and the lower sign 
corresponds to particles and antiparticles, respectively. The spin factor s = 1 (— 1) 
for 2;-spin up (down) is related to the helicity A by s = A sign(p^). Note that eq. (jHH)) 
is the dispersion relation in a general Lorentz frame, in contrast to the one derived 
in ref. 0. For the group velocity of the WKB wave-packet one obtains 

^- (l ± s--^ 
^0 V 2 EqEqz 

The semiclassical force acting on the particles, 

^ _ {m^y ^ ^ {m^ej ^ J'm^m^y ^^^^ 
2Eo 2EoEoz AEqEqz 

results from the canonical equations of motion. It was the main result of ref. ^3] 
that the expressions for the dispersion relation the group velocity (fHIIjl . and 
the semiclassical force agree with the full Schwinger-Keldysh result j23i . 

In the semiclassical approximation the evolution of the particle distributions /, is 
described by a set of classical Boltzmann equations. We assume a planar wall moving 
with constant velocity f„. Hence, in the rest frame of the wall, the distributions fi 
only depend on z, and p = |p|, due to the translational invariance parallel to the 
wall. For each fluid of particle type i we have 



[V 



gA + Fzd,Jf, = C,[f], (41) 

without any explicit time dependence, as we are looking for a stationary solution. 
The Cj are the collision terms describing the change of the phase-space density by 
particle interactions that drive the system back to equilibrium. We introduce pertur- 
bations around the chemical and kinetic equilibrium with the fluid-type truncation 
in the rest frame of the wall [H] 

M^^P^^P) = ^PME.+.L)-,.] ± 1 + ^f^^'^P^^P) (42) 



where (3 = 1/T and 'y^ = 1/ ^yl — v'^, and plus (minus) refers to fermions (bosons). 
Here the chemical potentials Hi{z) model a local departure from the equilibrium 
particle density and the perturbations 6fi describe the movement of the parti- 
cles in response to the force. The latter do not contribute to the particle density, 
i.e. J (Pp 6fi = 0. To first order in derivatives the perturbations are CP-even and 
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equal for particles and antiparticles. But to second order they have CP-even and 
CP-odd parts, which we treat separately, i.e. 



(43) 



so that the perturbations to second order for particles differ from those for antipar- 
ticles. 

In order to compute the asymmetry in the left-handed quark density, we expand 
the Boltzmann equation in gradients. In the model under consideration, the most 
important particle species are top and bottom quarks, as well as the Higgs bosons. 
The other quark flavors and the leptons can be neglected thanks to their small 
Yukawa couplings. In a first step we assume baryon number conservation. We take 
into account VT-scatterings, the top Yukawa interaction, the strong sphalerons, the 
top helicity flips and Higgs number violation with rates Pvf, Py, Pss, Pm and P/j, 
respectively, where the latter two are only present in the broken phase. After the 
left-handed quark asymmetry is computed, the weak sphalerons, with the rate P^g, 
convert it into a baryon asymmetry. 

We follow the computation and notation presented in ref. ^3]. We weight 
the Boltzmann equations with 1 and p^/i^o, and perform the momentum aver- 
age. Accordingly "plasma velocities" appear in the following, which are defined 
as Ui = {{pz/ Eq)5 fi) . We end up with the transport equations for chemical po- 
tentials of left-handed SU(2) doublet tops ^t,2^ left-handed SU(2) doublet bottoms 
/i6,2, left-handed SU(2) singlet tops /ifc,2, Higgs bosons /ift,2, and the corresponding 
plasma velocities 



-3P,,[(1 + 9K,^t)f^t,2 + (1 + 9K,^b)fib,2 + (1 - 9K,^t)f^t^,2] 







^VwKi^blJ,'^2 + 3Mfe,2 
— 3Pj/(/ib,2 + /^tS2 + fJ'h,2) — 3Piy(/^b,2 — fJ't,2) 
3P,,[(1 + 9ifi,0/it,2 + (1 + 9Ki,6)/ife,2 + (1 - 9Ki,t)/iic,2] 







3v^Ki^tfJ''tc^2 + 3t^„-ft'2,t("^?)W,2 + 3mJc^2 

-3Ty{j^t,2 + l^b,2 + 2/itc,2 + 2jJh,2) - 6Pr„(yUf,2 + l^t<=,2) 

3P,,[(1 + 9Ki,t)l^t,2 + (1 + 9Ki,b)l^b,2 + (1 - 9K,,t)l^t^,2] 







4v„Ki^hlJ''h,2 + 4Mft,2 
3Py(/ii,2 + f^b,2 + 2Aifc,2 + 2flh,2) - 4Pft/i/i,2 



(44) 
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= St 



-3i^4,t/iic,2 + 3v„K5,tUt-,2 + 3^;w^6,^("^?)'M^^2 + 3r*°*Mtc 2 = St 

-4i^4,h/i;,2 + 4t;wi^5A<2 + M/,,2 = 0. (45) 

Here the second-order perturbations label the difference between particles and an- 
tiparticles, i.e. ^2 = ^^20 — ^20 and U2 = U20 — M2o- On the r.h.s., St denotes the 
source term of the top quark, 

St = -v„Ks{mXy + v^K^e[ml{m])' . (46) 

The source term of the bottom quark, which is suppressed by ml/rn^ ~ 10"'^, has 
been neglected. The Higgs bosons do not have a source term to second order in 
gradients. The various thermal averages Ki in eqs. and are defined 

similarly to ref. We include the position dependence of the Ki. The damping 
of Ui^2 can be approximated by the total interaction rate, r*°*. In the numerical 
evaluations we have included a term 3Tw{ut^2 — Ub^2) which affects results only at 
the few percent level. Contrary to the transport equations in ref. |il4j we have 
doubled the degrees of freedom of the Higgs bosons to account for the second Higgs 
doublet in the model. 

Using baryon number conservation, the chemical potential of left-handed quarks 
can be expressed in terms of the solutions of the transport equations fit ^2, fJ'b,2 and 

/^t'=,2, 

1, , 
^^BL — /^gi,2 + /^g2,2 + 2 '>/^*.2 + fJ'b,2) 

= ^(1 + m^t)fit,2 + ^(1 + m,b)fib,2 - 2Kl,t/ii.,2. (47) 

Now, in a second step, the weak sphalerons convert the left-handed quark number 
into a baryon asymmetry. 



6 The baryon asymmetry 

The baryon asymmetry is obtained by 

r^s is the weak sphaleron rate, which is only present in the symmetric phase, and 
= 106.75 is the effective number of degrees of freedom in the plasma. The 
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Figure 6: The solid line represents 77^ as a function of the wall velocity for nih = 125 
GeV, tuh = 350 GeV, fij = 10000 GeV^ and = 0.2. This parameter setting 
determines = 4.5/T and ^ = 1.6. The dashed line would be the asymmetry when 
we substitute Eqz — > Eq in the dispersion relations. The dotted curve corresponds 
to the case where the ly-scatterings are in equilibrium. 

exponent p = 45r^s/(4fw) accounts for the relaxation of the baryon number in case 
of a slowly moving wall. 

In our evaluation we use the values F^^ = 1.0 x 10^® for the weak sphaleron rate 
F^, = 4.9xl0"^rforthe strong sphaleron rate ^,Ty = 4.2x lO^^T for the top 
Yukawa rate [Hlj, F^ = mf{z,T)/{63T) for the top helicity flip rate |54,, and = 
m^(z, T)/(50T) for the Higgs number violation rate [54]. Furthermore the total 
interaction rate can be expressed by the diffusion constant, F*°* = {DiKi i) / K^^i, 
where the quark diffusion constant is given hy Dq = 6/T [53] and the Higgs diffusion 
constant by Dh = 20/T 9 . The finite VT-scattering rate we approximate as Tw = 
-ptot rjj_^g bottom quark and the Higgs bosons are taken as massless. 

Fig. El displays the baryon asymmetry as a function of the wall velocity for one 
typical parameter set. The solid line indicates the solution when using the full set 
of transport equations (jl^ and If we resubstitute Eq^ — >• Eq in the dispersion 
relation the group velocity (jH^ and the semiclassical force (HO)), i.e. going 

back to these quantities as determined in ref . [HI , the resulting baryon asymmetry is 
substantially reduced (dashed line). This confirms the recent result that performing 
the boost back to a general Lorentz frame has a sizable effect and should not be 
neglected [Hj. In addition we have improved the transport equations by keeping 
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a finite VT-scattering rate. If tliese interactions were in equilibrium, rjB would be 
considerably overestimated for t>„ < 0.1 (dotted curve). We could also show that 
taking the Higgs bosons into account or not does not play a significant role. The 
same holds for the source terms proportional to the first-order perturbations fit,i 
and Ut^i, which we have neglected in the current paper. Their effect on the total 
baryon asymmetry is less than 10% in the model under consideration. 

Let us finally discuss the dependence of rjB on the Higgs masses. Fig. [7| shows 
contours of constant baryon asymmetry in the nih-mH plane, where we have fixed 
the parameters /i| = 10000 GeV^ and = 0.2. For each mass combination we deter- 
mine all quantities relevant for the phase transition, such as ^, tan Pt, Ogym and 
^brk to put them into the transport equations. There is only a mild dependence 
of the baryon asymmetry (cf. fig. El; so we consider only one wall velocity, = 0.1. 
In addition, the (^=l)-contour of fig. |21is also shown for orientation. As we increase 
TTiH, leaving rrih fixed, the asymmetry becomes larger. This behavior results from 
the mf ~ dependence of the top source term. Accordingly the baryon asymmetry 
becomes larger for a stronger phase transition. If we increase rrih, leaving the heavy 
Higgs mass fixed, rjB becomes smaller and reaches a minimum at rrih ~ 150-160 
GeV, similar to the behavior of L^. But in general there is only a minor depen- 
dence on the light Higgs mass. In this parameter setting it is possible to generate 
the observed baryon asymmetry for a heavy Higgs mass between 320 and 330 GeV 
and a light Higgs mass up to 160 GeV. Since rjB is more or less proportional to the 
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CP-violating phase 0, the measured value can also be explained for other values of 
the parameters if we adjust 0. Then the heavy Higgs mass should be somewhat 
larger. 

Comparing figs. El and we can use the baryon asymmetry to predict the EDMs. 
We see that for ^ ~ 1 and rrih = 115 GeV the neutron EDM is a factor of about 
2 above the experimental bound, which including the theoretical uncertainties is 
marginally tolerable. The electron EDM is a factor of about 5 below the experimen- 
tal bound. Moving along the = l)-contour to the largest Higgs mass of 190 GeV, 
we find \de\ ~ 0.03 x 10~^'^e cm and \dn\ ~ 0.7 x 10~^^e cm. Finally, taking ruh = 190 
GeV and mn = 400 GeV, we find |4| ~ 3 x lO'^^e cm and \dn\ ~ 0.09 x lO'^^e cm. 
So the experimental bound on the neutron EDM starts to cut into the parameter 
space of the model. Improving the bound by an order of magnitude would probe 
the larger part of the parameter space. The electron EDM is typically one to two 
orders of magnitude below the current bound. 

In this paper we focused on the case tan/3 = 1. As we discussed in the context 
of eq. ()28|) . larger values of tan/5 will lead to a smaller value of the baryon asym- 
metry, since the change in 6 is then mostly due to a change in 6i rather than 62. 
Extrapolating from the example of fig. El we estimate that for tan/5 > 10, successful 
baryogenesis should in any case be in conflict with the EDM bounds. It would be 
interesting to check this issue by direct evaluations. 

So there exists a wide range of realistic parameters where the computation of 77^ 
is under control, and which yields the observed baryon asymmetry. 

7 Conclusions 

We have studied electroweak baryogenesis in the 2HDM, focusing on the case of 
tan j3 = 1 and degenerate extra Higgs states. Evaluating the thermal Higgs potential 
in the one-loop approximation, we find a first-order phase transition, which is strong 
enough to avoid baryon number washout. This is achieved by the loop-contributions 
of the extra Higgs states, provided they are sufficiently strongly coupled. Taking 
/i3 = 10000 GeV^, this happens for a heavy Higgs mass thh > 300 GeV. The mass 
of the light, SM-hke Higgs, rrih, can be up to 200 GeV, or even larger. The Higgs 
potential allows the introduction of a single CP-violating phase, which has only a 
minor impact on the strength of the phase transition. In our example, if ttlh reaches 
about 500 GeV, the phase transition becomes very strong, while the perturbative 
description starts to break down. These findings are in agreement with those of 
ref. PI]. 

We have computed the properties of the phase boundary. The walls are typically 
thick, but the width decreases with larger tjih from L„ ~ 15T~^ to about 2T~^ We 
also compute the profile of the relative complex phase between the two Higgs vevs, 
which changes by an amount A6 between the broken and the symmetric phase. 

This phase shift leads to a CP-violating source term for the top quark, which 
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drives the generation of the baryon asymmetry. We compute the source term in the 
WKB approximation and solve the resulting transport equations, using the formal- 
ism of ref. 113]. We find that for typical parameter values the baryon asymmetry is 
in the range of the observed value. The explicit CP phase in the Higgs potential has 
to be taken between 10~^ and unity. For larger values of mn the baryon asymmetry 
increases, as the phase transition becomes stronger and the wall thinner. Our result 
differs from those of ref. [20] • There the baryon asymmetry in the 2HDM was com- 
puted using the method of reflection and transmission coefficients. In the regime of 
thick walls, this method is known not to give the leading contribution to the baryon 
asymmetry, which explains the different results. 

We have also computed the EDMs of the electron and neutron. Since there is 
only a single complex phase in the model, we can predict \de\ and \dn\ in terms of 
the baryon asymmetry and the Higgs masses. We find that > 10~^'^e cm. For 
the smallest allowed values of m/j and m/^, \dn\ can slightly exceed the experimental 
bound. Improving the neutron EDM sensitivity by an order of magnitude would 
test a substantial part of the parameter space of the model. The electron EDM 
is typically one to two orders of magnitude below the bound. These values are 
for tan/3 = 1. Extrapolating our results suggests that for tan (3 > 10, the 2HDM 
cannot produce the observed baryon asymmetry without being in confiict with the 
EDM constraints. In any case, the 2HDM can explain the baryon asymmetry for a 
considerable range of the model parameters. 

It would be interesting to extend our investigations to cover the full parameter 
space, in particular the case tan/? > 1. Since for larger values of uih higher- 
order corrections to the effective potential become more and more important, these 
contributions should be studied in more detail, most reliably on the lattice. Our 
proposal is testable at the LHC in the sense that at least one Higgs state should 
be observed. Furthermore, CP violation may be detectable in top pair production 
jSSEIl- Stringent tests could be performed at a future e~^e~ linear collider, where 
for instance deviations in the Higgs self-coupling could be detected jSHl- 
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